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SUMMARY 
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Explicit guidance equations are developed for steering multistage launch vehicles to 
injection conditions that satisfy typical lunar impact requirements. The calculus of vari- 
ations is used to derive the required boundary conditions on the thrust direction at staging 
points. An approximate variational yaw steering law is also derived in order to minimize 
the payload loss in removing yaw velocity errors. The equations are programed for 
steering the Atlas-Centaur vehicle on the direct-ascent Surveyor mission. Results are 
presented that demonstrate targeting accuracy and payload capability (compared with op- 
timized open- loop trajectories). ///, 



INTRODUCTION 

An important factor in determining the success of a mission is the ability of the 
guidance laws to steer the vehicle during its propulsion phases so that vehicle constraints 
are not violated and the required burnout conditions are achieved. The guidance laws 
must be capable of steering the vehicle to its desired target over a wide range of flight 
conditions, called dispersions, which are caused mainly by deviations from nominal ve- 
hicle performance. In addition to the requirement for targeting accuracy, guidance equa- 
tions should deliver near-maximum payload capability under nominal and dispersed flight 
conditions. 

Many guidance laws have been used and discussed in the literature. Reference 1 
contains a description of some of the most widely used formulations. With perturbation 
guidance schemes (e. g. , ref. 2), the thrust direction is calculated from empirical func- 
tions of the state conditions (position, velocity, and acceleration). These functions, and 
the associated constants, are obtained from curve fits of a nominal and a volume of dis- 
persed trajectories. The trajectories are flown with an optimized steering program 
(e. g. , calculus of variations) and iterated to the required terminal conditions (a process 
known as preflight targeting). Required velocity equations have also been used, wherein 



the thrust is pointed nearly parallel to the difference between the required and present ve- 
locity vectors. This method still requires empirical biasing functions to satisfy position 
constraints along the flight and to shape the trajectory to an optimum one. Again, pre- 
flight targeting is required to obtain the empirical functions and constants. 

The explicit guidance equations used in this report differ from the perturbation and 
required velocity equations in that the steering laws are based on an approximate closed- 
form solution to the equations of motion, solved continuously during flight. Several ad- 
vantages are inherent with this formulation: 

(1) The equations are not tied to any nominal trajectory so that no retargeting is re- 
quired when last-minute changes in vehicle or target parameters are introduced. 

(2) The equations have great flexibility and can be used without change for a variety 
of missions and vehicles. 

(3) The pitch and yaw steering laws closely approximate the form prescribed by the 
calculus of variations for maximum payload. 

The pitch steering law used in this report was first introduced by MacPherson 
(ref. 3), whose results were later expanded by Cherry (ref. 4). In reference 3, explicit 
equations are developed for steering a single stage, with constant thrust and specific im- 
pulse, from any initial position and velocity to a prescribed orbit, expressed in terms of 
radius, velocity, and flight-path angle. This pitch steering law allows a closed-form so- 
lution for radius and radial velocity (as functions of time). The required thrust direction 
is calculated directly as a function of the present and desired terminal state and the time 
to cutoff (determined by the required velocity). An approximation is developed to esti- 
mate the time to cutoff. Reference 4 uses the same pitch steering law, but extends the 
results to include other possible burnout requirements. An explicit solution is also de- 
veloped for engines that can be throttled. 

In this report, equations are developed that allow the explicit guidance laws to steer 
an arbitrary number of stages, each operating at a constant thrust and specific impulse. 
The calculus of variations is used to determine the required boundary conditions at 
staging points. In addition, the variational method is used to develop an approximate op- 
timum yaw steering law. A pseudotarget concept is introduced to compensate for Earth 
oblateness and other perturbing forces. 

The explicit equations are applied to the three-stage direct-ascent Atlas- Centaur- 
Surveyor mission, where the final two stages are flown with guidance- generated (closed 
loop) steering commands. The booster stage is flown through the atmosphere with a pre- 
specified (open loop) steering program. The powered flight trajectory places the vehicle 
on a coast ellipse that results in lunar impact at the desired location with a prespecified 
time of flight. In addition, the perigee radius of the transfer ellipse is specified. Equa- 
tions are developed that approximate these burnout requirements in terms of radius and 
radial velocity and thus allow the use of the closed-form solution obtained in reference 3. 
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The Surveyor mission has been chosen for study in this report because it represents 
an extremely difficult guidance problem. Among the mission and guidance requirements 
are 

(1) Closed- loop steering for two stages (sustainer and Centaur) 

(2) A velocity increment of approximately 27 000 feet per second added while steering 

closed loop 

(3) Variable launch azimuth and injection true anomaly, and the associated wide range 

of trajectory profiles 

(4) Accurate targeting to achieve lunar impact at the desired landing site with a pre- 

specified time of flight 

The explicit equations have been programed on an IBM 7094 computer along with two 
different vehicle and flight-environment simulations. A simplified program is used to 
obtain targeting and payload results, wherein ideal gravity and propulsion models are as- 
sumed. A detailed flight- simulation program is used to obtain dispersion results. The 
pseudotarget concept is used with this program to compensate for Earth oblateness and 
other perturbing forces. Targeting errors are expressed in terms of errors in injection 
perigee altitude, yaw velocity, and midcourse-correction requirements. 


ANALYSIS 

Explicit guidance equations are developed for steering the final two stages of a three- 
stage Atlas- Centaur launch vehicle to a set of injection conditions that result in a lunar 
impact at a prespecified landing site. The time of flight is also fixed in order to allow 
the Goldstone tracking station to view the lunar arrival. Because of atmospheric con- 
straints, the booster stage is flown with an open-loop pitch program designed to minimize 
aerodynamic heating and loads. The upper stages (sustainer and Centaur) operate under 
near-vacuum conditions, and both of these stages are flown closed loop. 

The following simplifying assumptions are used in the derivation of the guidance 
laws: 

(1) Spherical Earth, with an inverse- square force field 

(2) No atmospheric effects on upper stages 

(3) Constant-thrust specific- impulse operation for sustainer and Centaur stages 

(4) No perturbing bodies (in particular, a massless moon is assumed) 

The effect of these assumptions will be considered later in the analysis, and the equa- 
tions will be modified where required. 

One stage is assumed in the initial derivation; the required continuity equations at 
staging points will be considered later in the analysis. 
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Pitch Steering Law 


The steering laws are derived in a rotating radial, normal, circumferential coordi- 
nate system, defined by 


r = 


h = 


r x v ^ 


r x v 


9 = h x r J 


( 1 ) 


(All symbols are defined in appendix A. ) The assumptions stated earlier are used to 
write the differential equation for the radius 

't + -tL - co^r = a(f • r) (2) 
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where 


a(t) = 
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T 




T = 


a(0) 
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( 3 ) 


and ? is the unit thrust direction. In equation (3) and in the development that follows, 
t = 0 refers to the present time along the trajectory, rather than any fixed value of time. 
The guidance law used for the radial component of thrust is (ref. 3) 
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f • r = A + Bt +■ 


( 4 ) 


which is combined with equation (2) to give 
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¥ = (A + Bt)a (5) 

Equation (5) can be integrated twice to give 


r(T) = r(0) + b Q A + b^B 

(6a) 

r(T) = r(0) + r(0)T + c q A + CjB 

(6b) 


where 
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(7a) 


(7b) 


(7c) 


(7 d) 


and s is a dummy variable of integration. If the value of T (again, referenced to the 
present time) is known, equation (6) can be solved for values of A and B that will give 
the desired values of r(T) and r(T) (provided that these values are consistent with vehi- 
cle capability). For the Surveyor mission, none of the variables T, r(T), or i-(T) is 
explicitly specified. In order to make use of the closed-form solution, these variables 
must be specified in terms of the required Surveyor boundary conditions. 


Surveyor Boundary Conditions 

The (prespecified) translunar flight time essentially depends only on the injection 
energy, so that the injection energy must also be fixed. The shape of the transfer ellipse 
is then completely determined by the perigee radius. Maximum payload capability is ob- 
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tained by injecting at low altitudes (and perigee altitudes), which causes payload heating 
constraints to be violated. For this reason, the perigee radius has been fixed for Survey- 
or flights. The values of energy and perigee radius determine the eccentricity and semi- 
latus rectum of the translunar ellipse: 


e 




(8a) 


p = JL (e 2 - 1) 
2E 


(8b) 


If the required true anomaly t](T) were known, the values of r and r at cutoff could be 
calculated from 


r(T) = 2 (9a) 

1 + e cos ri (T) 

r(T) = ^ e sin 7] (T) (9b) 

For Surveyor trajectories, cutoff occurs when the required energy level has been 
achieved. The required true anomaly is determined by the Earth-moon geometry at im- 
pact and the length of the boost arc from launch to injection. Since the cutoff time and 
the length of the boost arc are not precisely known during flight, the explicit guidance 
equations require estimates of T and tj(T) in order to determine r(T), r(T), and the 
steering constants A and B (from eq. (5)). These estimates should be reasonably ac- 
curate early in flight and should converge rapidly near cutoff, so that accurate targeting 
can be achieved with no appreciable degradation in performance. The equations used to 
estimate T and r](T) will be developed later in the analysis. 

Yaw Steering Law 

Assuming an inverse- square gravity field (as stated earlier) gives the required in- 
jection condition on yaw velocity as 



(10) 
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Figure L - Lunar transfer geometry. 


where R is a unit vector directed 
m 

at the target, assumed fixed in space 
(see fig. 1). It is stated in equation 
(10) that the injection r,v plane 
contains the target vector. Differen- 
tiating equation (10) gives 

-(h-R m ) = - — (f-h) ( 0'H m ) 

dt m Vq 

( 11 ) 

where 

V 0 = 7- 9 

Equations (10) and (11) can be used 
to restate the requirement on yaw 
velocity: 



— (f- m • R m )dt = 0 
V 6 


( 12 ) 


A component of thrust must be supplied in the yaw direction to remove the initial yaw ve- 
locity. Moreover, the yaw steering program must be optimized so that the payload loss 
due to yaw steering is minimized. The appropriate yaw steering law is derived in appen- 
dix B by using an approximate calculus- of- variations solution: 


f- 



(13) 


The gain constant K can be calculated by combining equations (12) and (13): 



(14) 
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The integral in equation (14) cannot be explicitly determined during flight, so that an ap- 
proximation is required to determine the gain constant K. 


Yaw Integral Approximation 


In order to evaluate the integral in equation (14), a linear variation is assumed for 


<«' R m>/V 


e • r 


m 


d l + d 2 4 


(15a) 


where 



(15b) 


(15c) 


Equation (7) is now used to evaluate the yaw integral 


which gives 


and 
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(h- L) 

\ m / 0 


( d l b 0 + 2d l d 2 b l + d 2 b 2) 


f • h = 


( h -R m )(0-R m ) 


v e ( d i b o + 2d i d 2 b i + d 2 b 2 ) 


(16) 


(17) 


(18) 
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Estimate of Required True Anomaly 


The differential equation for polar angle can be written 


e = - (f ■ e) - — 


( 19 ) 


which is integrated twice to give 


A6 = e (t) - e ( 0 ) = 


Mil 


(f • 0)ds 


dt - 2 


m 


— ds ldt (20) 


Since the variables f • 9 and h are not known explicitly as functions of time, several 
approximations are required in order to evaluate the integrals in equation (20). The 
thrust vector is directed close to the local horizontal in order to increase energy and an- 
gular momentum. Thus, the components of thrust along the radial and normal directions 
are relatively small so that the approximation 


f • e = Vi - (f- r) 2 - (f- h) 2 ~ l - r > 2 - 


( 21 ) 


may be used. In addition, f • r and f • h (from eqs. (4) and (18)) can be assumed linear 
in time: 


f • r ~ f + ft 
r r 


(22a) 
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co ri 


f r = A +| 


(22b) 
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f r = B + 


(22c) 


9 


and 
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< h ' R m> d l 


£ _ A 

( d l b 0 + 2d l d 2 b l + d 2 b 2) 


f - — f 
h " H h 

d l 


Using these approximations gives f • 9 in the form 


f • e = fg + tgt + fgt^ 


where 


f 2 f 2 

0 2 2 


1 9"-(«r i r +1 h^ 


f 0 = - 


;2 j2 

f r + f h 


(23a) 

(23b) 

(23c) 

(24) 

(25a) 

(25b) 

(25c) 


The radius varies only a few percent during the flight, so that r may be assumed con- 
stant in the first integral in equation (20): 

r ~ 7 = r( ° ) + r(T ) (26) 
2 

o 

The expression 2hr/r is the Coriolis acceleration, which is small compared with the 
thrust acceleration. A linear variation is therefore adequate for this term: 

— ~ dg + d 4 t (27) 
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where, as previously, 



(28a) 


(28b) 


The preceding approximations and equation (7) are used in the integration of equation (20) 
to give 


A0 



+ f 0 c l + f 0 c 2 


) - 


d 3 T 


d 4 T 3 

3 


(29) 


Equation (29) leads to an estimate of the required true anomaly. From figure 1 

r](T) = T] m + A6 - d (30) 

where 77 m is the target vector true anomaly. 


Estimate of Cutoff Time 


The required angular momentum at cutoff is given by 


h(T) = VMP 


and the time derivative of angular momentum is 


h = ar(f- 9) 


The required change in angular momentum is 


Ah = h(T) - h(0) 


-c 


ar(f • 0 )dt 


(31) 


(32) 


(33) 
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The approximations for f • 6 and r from equations (24) and (26) are used to give equa- 
tion (33) in the form 


Ah = r(f 0 b o + f 0 b i +? 0 b 2 ) 


(34) 


Equation (34) can be expressed in terms of ideal Av by expanding b. Using equation (7) 
yields 


Ah - r 


(f 0 + f 0 T + f 0 T 2 )AV-V e T(f 0 +f 0 T) - 




(35) 


or 


Av = 


^ + v e T ^ + V) + — 


h + V + V 2 


(36) 


The terms containing T in equation (36) are small compared with Ah/r so that the up- 
dated estimate of T from the previous computer cycle can be used in evaluating Av. 
The new estimate of T is obtained by combining equations (7a) and (36): 

Av = -v e log ^1 - — j (37a) 

/ -Av/v \ 

T = t( 1 - e e ) (37b) 


Extension to Two Stages 

A A 

The calculus of variations shows (ref. 5) that the thrust direction and rate f and t 
should be continuous across staging for optimum performance. Since the acceleration 
will be discontinuous at staging points, A and B must also be discontinuous in order to 

A -A 

satisfy the continuity of f and f. Differentiating equation (4) results in 
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— (f- r) = B 
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- 2 o)(f- e) 


(38) 


so that the required continuity equations are 
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and the discontinuities in A and B are 


l 2 ' A 1 " °1 A 1 


(jL - a) 2 r\ 

1 1 

l 2 l 

aj(T^) a 2 (0) 


(40a) 


^r 


)(i iV 

Lj _ mV 


1 1 

/ g \ v e,l v e,2/ 

L\ r 2 ) J 

s 

a-j(Tj) ^(O) 


(40b) 


The additional subscripts 1 and 2 introduced in equations (39) and (40) refer to the 
sustainer and Centaur stages, and the subscript s refers to conditions at staging. The 
time reference point for the Centaur stage is set at Centaur startup (during sustainer 
phase) . 

For two stages, equation (6) becomes 
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r(T) - r(0) + (bg ^ + bp 2 )-Aq + ^1 1 + ^1 2 + b 0 2^'l^ B l + ^0 2 + ^*1 2 AB (41a) 

r (T) = r(0) + r(0)(T 1 + T 2 ) + (Cg, i + c Q 2 + b^ j^^l 


+ < c l, 1 + b l, 1 T 2 + c 0, 2 T 1 + c l, 2^ B 1 + c 0, 2 AA + c l, 2 AB < 41b > 

In order to determine AA and AB from equation (40), estimates of position and ve- 
locity at sustainer cutoff are required. These estimates are obtained by using equations 
(6) and (34). 


)) + b Q jAj +b 1>1 B 1 

(42a) 

’(0)T 1 + c Q , jAj + c^ jBj 

(42b) 

(f 0 , l b 0, 1 + V l b l, 1 + * 6 , l b 2, 1> 

(42c) 

h(T-i) 

(TJ = ‘ 

1 r (T-i) 

(42d) 


Boundary Condition on Yaw Thrust Direction 


The calculus of variations requires the yaw thrust direction to be continuous across 
staging, which implies that K in equation (13) is continuous. However, the derivative 
of f • h is not continuous, since the derivative of h involves acceleration (and is conse- 
quently discontinuous) . Thus, different values of d^ and d 2 are required for each 
stage. 

Extension of Estimates to Two Stages 

Yaw integral . - For two stages, equation (16) becomes 



Boost arc . - The estimate of A 9 for each stage is 


14 


49 f = 


— ] T. + — (f 


! c n a + ^ n c i n + /iC 
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and the total A0 is 


A0 = A0j + A#2 (45) 

Cutoff time . - The sustainer stage cutoff time is assumed fixed, equal to its nominal 
value. Equation (42c) is used to estimate the angular momentum at Centaur startup, and 
equation (36) supplies the estimate of second stage Av, with Ah = h(T) - h(T^). Equa- 
tion (37) is then used to estimate the Centaur stage cutoff time. 

The staging logic presented in the preceding section is, of course, not limited to 
two stages. Equations (38) to (45) are derived for an arbitrary number of stages (N) in 
appendix C. 


Cutoff Logic 

As cutoff energy is approached, the equations for calculating A, B, T, and K be- 
come indeterminate. If these calculations are continued up to cutoff, the values begin to 
oscillate and eventually diverge. For this reason, it is necessary to terminate the major 
calculation loop prior to cutoff and to steer the vehicle for the rest of the flight by using 
the last calculated values in equations (4) and (18). It is important, however, not to ter- 
minate these calculations too far from cutoff or targeting accuracy may be degraded be- 
cause of errors in the last estimates of A and B. For Atlas- Centaur- Surveyor simula- 
tions, a value of "energy-to-go” corresponding approximately to 10 seconds prior to 
cutoff has been used with good results to terminate the major calculation loop. 

Although the last estimates of A and B allow accurate targeting, the last estimate 
of cutoff time is not accurate enough to use for actual termination of powered flight. A 
more accurate method has been developed for this purpose in which energy is assumed 
to be a parabolic function of time. The coefficients of the energy-time polynomial are 
obtained by using the last three calculated time and energy values prior to cutoff. The 
cutoff time is determined by solving the quadratic. 


EXTENSION TO REAL WORLD PROBLEM 

The equations in the preceding sections have been derived by using several simplify- 
ing propulsion and gravity assumptions. Adjustments for these assumptions must be 
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made in two areas - navigational equations and steering equations. Both these areas will 
be considered in the following sections. 


Propulsion Effects 

The thrust and specific impulse of a real vehicle are not constant throughout flight, 
but the variations are usually small enough so that equation (3) can be used without de- 
grading performance or accuracy. The measured acceleration and nominal specific im- 
pulse (for each stage) are used as inputs to the equations to estimate the future accelera- 
tion history. No in-flight measurement of specific impulse is assumed, because the 
noise in such measurements generally exceeds the expected (statistical) deviation from 
the nominal value. Simulations have shown that targeting accuracy is not degraded when 
large 3a dispersions in specific impulse are introduced. 


Atmospheric Effects 

Atmospheric forces are large during booster phase, but this stage is flown open- loop 
so that the guidance equations are not affected. All contact forces are detected by the 
accelerometers so that navigational accuracy is also unaffected. Since atmospheric 
forces have greatly diminished at guidance closure, no special logic is needed in the 
steering equations to compensate for these effects. The guidance equations simply treat 
all measured forces as if they were a result of thrust. 


Gravity Effects 

In the preceding analysis, a spherical Earth and no perturbing bodies were assumed. 
The real Earth, however, is oblate and does not have a simple inverse- square gravity 
field. In addition, the translunar trajectory is perturbed by the moon, the sun, and the 
planets. Since the guidance accelerometers do not detect gravitational forces, a gravity 
model must be prestored on the guidance computer and used in the navigational equations. 

One possibility is to store an oblate Earth gravity model, as well as ephmerides for 
the moon, the sun, and the planets, and thus determine the trajectory with great preci- 
sion; however, this procedure is undesirable for several reasons. First, the free-flight 
equations of motion cannot be solved analytically when perturbing forces are present, so 
that the required injection conditions cannot be precisely determined. Second, storage 
limitations on most flight computers would preclude the possibility of storing such a 
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large amount of gravity data. 

An alternate method that has been used with good results is the pseudotarget method, 
wherein a spherical gravity model is assumed for Earth and no perturbing forces are con- 
sidered. An iteration is performed in the targeting process to determine the direction of 
the pseudotarget vector, as well as the pseudoenergy and perigee radius, which satisfy 
the real world injection and lunar impact requirements under nominal flight conditions. 
When dispersions are introduced, the trajectory profile varies somewhat, causing the 
compensated perturbing forces to have a slightly different effect on the trajectory. Tar- 
geting accuracy is somewhat degraded under these conditions; however, the magnitudes 
of the perturbing forces are small, so that their variation on perturbed flights is a minor 
effect, and targeting accuracy is not degraded significantly. 


SURVEYOR MISSION 


The direct-ascent Surveyor launch windows are characterized by the systematic 
variation of launch azimuth and injection true anomaly with time, as shown in figures 2 
and 3. The variation of launch azimuth and true anomaly through any launch window is 
caused by the rotation of the Earth, while the particular curves to be followed in figures 
2 and 3 are dictated mainly by the declination of the moon at arrival, as shown in the fig- 
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Launch azimuth, deg 


Figure 2. - Launch azimuth as function of time of launch. 


Figure 3. - True anomaly as function of launch azimuth. 
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Figure 4. - Variation of burnout weight capability with injection Figure 5. - Variation of burnout 

true anomaly. Launch azimuth, 108 degrees. weight capability with launch 

azimuth for fixed true anomaly. 


ures. The choice of zero time in the figures is arbitrary. 

The duration of each launch window is determined by range safety launch azimuth 
limits, as indicated in figures 2 and 3, and by minimum allowable payload weight. The 
variation of payload capability with true anomaly and launch azimuth is shown in figures 
4 and 5. The specification of a minimum payload weight limits the allowable true anom- 
aly range, as shown in figure 4. The explicit guidance equations must be capable of 
steering the vehicle within the allowable range of launch azimuth and true anomaly. 


VEHICLE DEFINITION 

The Atlas-Centaur launch vehicle is sometimes referred to as a two and one-half 
stage vehicle because of the use of the sustainer engine during booster phase as well as 
during sustainer solo phase. In addition, the booster and sustainer engines share pro- 
pellant tanks. 

The vehicle lifts off with the booster and sustainer engines operating. At a predeter- 
mined acceleration level, the booster engines are shut down and jettisoned, while the 
propellant tanks are retained for use by the sustainer engine. Sustainer engine cutoff oc- 
curs at propellant depletion, at which time the sustainer engine and propellant tanks are 
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TABLE I. - VEHICLE AND FLIGHT ENVIRONMENT ASSUMPTIONS 


Performance 

parameter 

Simplified simulation 

Detailed simulation 

Booster 

Sustainer 

Centaur 

Booster 

Sustainer 

Centaur 

Thrust 

F sf +A ex AP 

Constant 

Constant 

Detailed propulsion 
model 

Constant 

Weight flow 

Constant 

Constant 

Constant 

Constant 

Atmospheric 

forces 

Yes 

No 

No 

Yes 

Yes 

No 

Gravity 

forces 

Spherical Earth, no 
perturbing bodies 

Oblate Earth model, plus 
perturbing forces due to 
sun, moon, and Jupiter 

Steering 

Zero 
angle of 
attack 

Closed loop 

Series 

of 

constant 

pitch 

rates 

Closed loop 


jettisoned. After a short delay, the Centaur engines ignite and burn until mission energy 
is achieved, at which time the Centaur engines are shut down. 

Two different vehicle and flight environment simulations are used in the guidance 
equation simulation studies. The assumptions used with the two simulations are outlined 
in table I. 

The detailed simulation uses a detailed propulsion model for the booster and sus- 
tainer engines. With this model, booster and sustainer thrust and weight flow are calcu- 
lated as reference values plus linear perturbation terms that depend on atmospheric con- 
ditions, pump inlet conditions, acceleration, etc. The Centaur engine is assumed to op- 
erate at constant thrust and weight flow. A detailed atmosphere model is also used with 
this program, wherein aerodynamic lift and drag are considered and gimbaling of the 
booster and sustainer engines is used to cancel aerodynamic moments. The gravitational 
model includes a detailed oblate Earth and perturbing forces due to the sun, moon, and 
Jupiter in order to determine a precision translunar trajectory. 

The simplified program uses only an exit-area - pressure correction term for 
booster thrust, while sustainer and Centaur thrust are assumed constant. All weight 
flows are assumed constant. The numerical values for this program are obtained from 
a "best fit" of the detailed trajectory. 

With the simplified program, the only atmospheric force considered is axial drag 
during boost phase. The Earth is assumed spherical, with an inverse- square force field, 
and no perturbing bodies are considered. 
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(a) Enter guidance program. (W Navigation and booster. 










(If sustainer phase 
exceeds expected 
duration, t must 
be made positive. ) 
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Figure 6. - Explicit guidance equations flow charts. 






















The simulation of the explicit guidance equations used with both the simplified and 
detailed programs is outlined in figure 6. The flight computer navigational equations 
are not simulated; the real and guidance trajectories are both integrated using a fourth- 
order Runge-Kutta integration scheme. Flight computer navigational errors can be 
assumed small and are ignored here. All contact forces are used in the guidance- 
integrated trajectory, but a spherical Earth gravity model is used, consistent with the 
pseudotarget concept. Vehicle and computer time lags and response characteristics are 
not simulated. 


RESULTS 

The targeting accuracy and payload degradation of the explicit equations are demon- 
strated by using the simplified program. No targeting is required to compensate for 
gravity force with this program, since an ideal gravity model is used in the trajectory 
integration. The results obtained with this program are presented in table n. 

The trajectories presented in table II cover the range of possible launch azimuth and 
true anomaly, as defined by figures 2 to 5. Targeting accuracy is evaluated in terms of 
errors in injection true anomaly, yaw velocity, and perigee radius. These errors are 
related to midcourse velocity requirements by approximate formulas derived in appen- 
dix D. Payload degradation is determined by comparison of guided trajectories with tra- 
jectories flown by using a calculus- of- variations steering program. 


TABLE H. - TARGETING AND PAYLOAD RESULTS 


Launch 

azimuth, 

deg 

True 

anomaly, 

deg 

Error in 
yaw velocity, 
ft/sec 

(a) 

Error in 
perigee 
radius, 
ft 
(a) 

Error in 
true anomaly, 
deg 

(a) 

Payload 

loss, 

lb 

(b) 

Midcourse 

correction 

requirements, 

ft/sec 

92 

15.4 

-0. 035 

-146 

0. 00104 

4.9 

0. 143 

112 

12. 3 

-.058 

-125 

. 00087 

4. 3 

. 110 

104 

9.9 

-.029 

-104 

.00053 

4.3 

.073 

99 

7.3 

-.021 

-90 

-. 00002 

4.2 

.007 

96 

4.6 

-.013 

-91 

-. 00005 

4.2 

.009 

94 

2.0 

-.005 

-89 

-. 00004 

4. 5 

.006 

92 

-0.8 

-.002 

-89 

-. 00004 

4.9 

.006 

91 

-3. 5 

.003 

-89 

-. 00096 

5. 4 

. 131 

95 

-6.4 

0 

-104 

-. 00048 

7.4 

.066 


a ActuaI value minus desired value. 

^Explicit guidance equations minus calculus of variations. 
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TABLE HI. - DISPERSION RESULTS (DETAILED PROGRAM) 


Dispersion 

Magnitude 

Error in 
perigee 
radius, 
ft 
(a) 

Change in 
weight, 
lb 

(a) 

Midcourse correction 
requirement, ft/sec 

Miss only 
(b) 

Miss + time 
(c) 

Booster specific impulse 

-3. 63 sec 

-6 

-81. 4 

0. 197 

0. 550 

Sustainer specific impulse 

-3. 07 sec 

23 

-51. 3 

.088 

. 148 

Booster staging 

-0. 08 g 

-26 

3. 1 

.931 

2. 830 

Booster pitch program 

5 percent 

-33 

-18.2 

.062 

. 370 

Booster pitch program 

- 5 percent 

21 

-15.0 

1.008 

2. 820 

Centaur thrust 

424 lb 

35 

-13. 2 

. 288 

1. 188 

Centaur specific impulse 

-4. 95 sec 

67 

-135. 0 

.206 

. 302 

Centaur expendables 

403 lb 

-15 

37. 8 

. 655 

1. 781 

Launch azimuth 

2 deg 

-30 

-. 6 

1. 399 

5. 260 

Launch azimuth 

-2 deg 

27 

-2. 3 

1.257 

3. 862 

Root sum square 


96 

171. 5 

2. 108 

6. 770 


a Dispersed minus nominal. 

^Midcourse AV to correct miss, allowing optimum time of flight. 
c Midcourse AV to correct miss and time of flight errors. 


In order to demonstrate the validity of the pseudotarget concept, the detailed simu- 
lation program is used to target a trajectory under nominal conditions, and various dis- 
persed trajectories are flown. These results are presented in table III. 

The dispersions used in table III are statistical 3a values. Although many other 
dispersions are possible, those presented in table III have the largest effect on injection 
conditions and midcourse velocity requirements, and the root- sum- square values ob- 
tained represent over 90 percent of the total obtained by using the complete dispersion 
list. 

The results in table III show that perigee control is excellent in spite of the approxi- 
mate gravity model used. The weight deviations are the result of performance losses 
associated with the dispersions, rather than degradation by the equations. The root- sum- 
square weight deviation represents the flight performance reserve required to compen- 
sate for 3a performance dispersions. 

Midcourse velocity requirements are calculated by using sensitivity coefficients, 
obtained by integrating the nominal trajectory and the adjoint equations to the moon on 
the detailed N-body program. Again, the root- sum- square values for miss only and 
miss plus time of flight indicate good performance. 

Typical pitch and yaw steering profiles generated by using the explicit equations 
with the simplified trajectory program are presented in figures 7 and 8. The optimum 
calculus- of-variations profiles are also presented for comparison. The discontinuities 
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observed at sustainer cutoff are caused by errors in estimating sustainer cutoff conditions, 
and by the time lag before Centaur startup (ignored in the equations). 


CONCLUDING REMARKS 

Explicit guidance equations have been developed for steering a three-stage Atlas- 
Centaur boost vehicle to injection conditions that satisfy the Surveyor mission require- 
ments. The targeting and dispersion results presented in tables II and III demonstrate 
the precision with which the explicit equations satisfy these requirements. 

Before the equations can be implemented in a real vehicle, several additional study 
areas are required: 

1. The explicit equations require more computer storage than the semiempirical 
velocity equations presently in use for the Surveyor mission. Preliminary estimates in- 
dicate that the explicit equations probably would not fit into the present Centaur computer, 
but would fit into several other existing and proposed flight computers. 

2. Because of the large number of calculations required, cycle time could become 
excessively long, dependent, of course, on the computer speed. If this occurs, major 
and minor cycles could be introduced in which the latest steering constants are used to 
supply steering signals in the minor cycle (steering block, fig. 6(h)), while the steering 
constants are updated less frequently in the major cycle. 

3. Some study is required to define (and, if necessary, to minimize) the coupling be- 
tween the steering equations and the vehicle due to computer and autopilot time lags and 
response characteristics. These effects are related to the computer cycle time and 
should not present a problem because the minor cycle time can be made arbitrarily short. 

In conclusion, the explicit equations presented appear to be an excellent choice for 
use with the direct-ascent Surveyor mission. The targeting and performance of the equa- 
tions are excellent, and no severe problem areas appear to exist. 

Because of the inherent accuracy and flexibility of the equations, their use can also 
be recommended for other lunar and planetary missions, both direct-ascent and parking 
orbit, as well as other less complex missions. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 6, 1965. 
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APPENDIX A 


SYMBOLS 


A 

steering constant (eq. (5)), 

K 

yaw steering gain constant, 

A ex 

nondimensional 

2 

exit area, ft 

M 

ft/sec 

target miss distance, ft 

a 

o 

acceleration, ft/sec 

m 

mass, slugs 

B 

steering constant (eq. (5)), 
sec - * 

P 

2 

pressure, lb/ft 


P 

semilatus rectum, ft 

b, c 

steering integrals (eq. (7)) 

*m 

unit target vector, nondimen- 

d 

coefficient of linear polynomial 

sional 


(eqs. (15) and (28)) 

r 

radius, ft 

E 

2 2 

energy per unit mass, ft /sec 

r 

m 

target distance, ft 

E c 

energy-to-go at which steering 

r 

perigee radius, ft 


constant calculations are 
2 2 

terminated, ft /sec 

P 

r 

average radius, ft 

e 

eccentricity, nondimensional 

r,0,h 

unit radial, circumferential, 

F 

thrust, lb 


normal coordinate system, 
nondimensional 

f 

unit thrust direction, nondimen- 

T 

cutoff time, sec 


sional 

T 

time of flight from midcourse 

f x’ f x’ k 

x- components of unit accelera- 

m 

to target, sec 

G 

eration direction (x = r, 6, h) 

effective radial gravity acceler- 
2 

ation, ft/sec 

CQ 

-M -W 

time, sec 

estimate of sustainer cutoff 

g 

Earth surface gravity accelera- 
tion, ft/sec 2 

V 

time, sec 
velocity, ft/sec 

h 

angular momentum per unit 

v e 

jet velocity, ft/sec 


mass, ft 2 /sec 

v h 

yaw velocity, ft/sec 

*sp 

specific impulse, sec 

Av L 

velocity loss due to yaw steer- 

J 

functional to be minimized, 


ing, ft/sec 


lb-sec 2 /ft 

v e 

circumferential velocity, ft/sec 
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«,/3,y 

parameters used to calculate 

d 

desired 

6, A 

steering constants (fig. 6) 

f 

final 

e 

2 2 

energy-to-go to cutoff, ft /sec 

m 

target 

e h 

yaw error, nondimensional 

r> 

s 

staging point 

X 

M 

Lagrange multiplier, lb- sec /ft 

sjf 

sea level 

M 

central body (Earth) force con- 
o o 

T 

cutoff 


stant, ft /sec 





0 

present value 

V 

true anomaly, rad 





1 

sustainer phase 

e 

polar angle, rad 





2 

Centaur phase 

<p 

angle from present position to 




target (fig. 1), rad 

Superscripts: 

T 

specific stage time, sec 

C) 

derivative of l, 




to time 

O) 

angular velocity, rad/sec 





— ► 

vector 

Subscripts: 






unit vector 

c 

midcourse 




with respect 
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APPENDIX B 


DERIVATION OF YAW STEERING LAW 


The problem to be solved is to zero the initial yaw velocity error while minimizing 
the payload loss due to yaw steering. This payload loss is a result of a reduction in ef- 
fective acceleration in the pitch plane. The amount of this reduction is 

Aa = ajl - Vl - (f • h) 2 j ~ | (f • h) 2 (Bl) 

which gives a pitch plane velocity loss of 


Av 


L " 



^(f- h) 2 dt 
2 


(B2) 


The velocity loss must be regained by expending additional propellants. Using equation 
(37a) results in 


Am T = m T 


( -Avj/v \ m T 
(l - e L e ) - Avj 


Equation (12) is used to state the problem mathematically as minimizing 


Am T 


/* T - (f • h) 2 dt 

-eJo 2 


(B3) 


(B4) 


subject to the constraint 



(f- h)(0-R m )dt=(h. R m ) Q 


(B5) 


This problem can be solved by using the calculus of variations and the method of undeter- 
mined Lagrange multipliers (ref. 6). Specifically, 
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J = 



m T a 



(f- h) 2 


dt + X 



a 



(f- h)(0. R m )dt 


(B6) 


where X is an undetermined Lagrange multiplier. The payload loss is minimized when 
J is minimized and equation (B5) is satisfied. Taking the differential of J yields 


6 J = 



(f- h) +X 


(6 • R m ) 


6(f • h)dt 


(B7) 


where a, (6 • R m ), and have been assumed definite functions of time. To minimize J, 
it is required that 8J = 0. Since the variation of (f • h) is arbitrary, the term in brackets 
must be set equal to zero. 


f- h = K 


'0 • R. 


m i 


(B8) 


where 


K= - 



(B9) 


Equation (B8) is the yaw steering law to be used. The multiplier K is evaluated as in 
equation (14): 


K = 



(BIO) 
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APPENDIX C 


DERIVATION OF N-STAGE CONTINUITY EQUATIONS 


The discontinuities in A and B at each staging point are (following eq. (40)) 


AA i = A i + 1 ‘ A i ' B i T i = (— - uj2r ) (— “ . 

V 2 /s \ a i <T J ) a j +1<0 V 

j 


and 


iB rVrV 


u 2 
— — U) i 



v V • 1 

e,J e,j+l 



a j (T j) a j+l (0) / 


(Cl) 


(C2) 


and the steering constants for each stage are 


B 


£=i 


(C3) 


and 


i-i 

V A i + Z ( AA * +B * T *) 

f=i 


i-i / f-i 

A 1 + Z r A f + B 1 T £ + Z AB k, 

£=1 \ k=l 


(C4) 


The change in r for each stage is 


4 r r b o,jV b i,j B j 


/ i- 1 \ i_1 / 

= b o,i A x + ( b i,j +b o,iS T t) B i \ b o,i iA ( +b o,i T t 2 j 

\ £=1 / i =1 \ k=l 


AB k +b l J AB £ 


(C5) 
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4 


and the initial r for each stage is 


j-1 

i'j(O) = i-j (0) + ^ Ar f 
£-1 


j-1/ £-1 

2\ b l,£ +b 0,£ 2 T t 

‘ ' k=l 


£- 1 ' 


j-1 £-1/ k-1 

B l + 2 2( b 0,£ AA k +b 0,£ T k2 AB i 
£-lk-l\ i = i 


+ b l, £ AB k 


The change in r for each stage is 

ar r r i ( 0 ) V c o,iV c i,i B j 


(C6) 


/ 3- 1 \ 

j-1 / £-1 

j-1 

( C 0, j +T j 2 b 0, £y A l + 

c i,i + E c o,j T « +b x,f T j +b o,« T jE T k 

b i + 2 c 

\ £=1 / 

£=1 \ k-1 / 

£-1 


J 0,j a/i £ 


+ C 1J AB £ + 


£-1/ k-1 x 

I b 0,t T j aA k +b 0,« T k T j E aB i +B l,t T j aB k + c o,i T t aB k 

k=l\ i-1 / 


+ r 1 (0)T j 

(Cl) 


and the initial r for each stage is 

j-1 

r.(0) = rj(0) Ar f 

£-1 


/ ^ \ j" 1 

^ / *-1 \" 

( c 0, £ +T £ 2 b 0,k) A l + 2 

c l, £ + 2 ( c 0,£ T k +b l, k T £ +b 0, k T £ 2 T i ) 

\ k-1 / £=1 

k-1 \ i-1 / 


j-1 

E 

£=1 

j-1 £-1 

+ EE 

£=1 k-1 


B, 


k-1 

c 0, £ AA k +c l,£ AB k + 2 ( b 0, k T £ AA i +b 0,k T i T £ 2 AB m 

i-1 


i-1 


m=l 


+ b l,k T £ AB i + c 0,£ T k AB i 


j-1 

-1(0)2 t £ +ri (0) 
£-1 


(C8) 
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The final r and r are obtained from equations (C6) and (C8) by substituting N + 1 


Bi 


N £-1 / k-1 

+ IE E b 0,f AB i +b l,£ A3B k) 

£ = lk=l\ i=l 


/ N \ 

N / £-1 \ 

r(T)=r 1 (0) + ( £ b Q J A 1 + 

2 ( b l, £ + b 0, £ El T k) 

\f =1 / 

£=1 V k=l / 


N / £-1 \ N 

r < T > = I c 0,£ +T £E b 0,k A l + S 


£-1 


k-1 

b 0, k T £ S T i 


i=l 


k=l / £=1 

N £-1 

B 1 + IZ 

£=1 k=l 


: 1,£ + X (^0,£ T k +b l,k T £ 

k=l \ 


k-1 


c 0,£ AA k + c l,£ AB k + El ( b 0,k T £ AA i 

i=l 


i-1 


+ b 0,k T i T £ El AB m +b l,k T £ AB i + c 0,£ T k AB i J 

m=l 


N 


+ r t (0) El T £ + r l(°) 

£=1 


The estimates of staging conditions are obtained by using 

r(Tj) = r J+1 (0) 

r <Tj) ' r j +1 <°) 

|r,(0) + r, .(0)1 . 

h(Tj , . -.(o) ♦ Li r l.i (f 9> ,b 0> , + , ♦ v ib 2> j) 


h(T,) 

v 6 (T.) = i- 

6 ] r(Tp 

The yaw integral in equation (43) becomes 



N 




£ + 2d l,£ d 2,£ b l,£ +d 2,£ b 2,j 


for j: 
(C9) 


(CIO) 

(Cl la) 
(Cl lb) 

(Cllc) 

(Cl Id) 

(Cl 2) 
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^ * 


and the travel angle for each stage is 


vH* 

c +f c +f c Vd T 2 _ d M T * 
M2,jj a 3, £ a £ 3 

(C13) 

' 0,1 1 




Equations (36) and (37) are used as in the text to estimate the final stage cutoff time, with 


Ah = h(T) - h N (0) 

The cutoff times of the first N - 1 stages are assumed fixed and equal to their nominal 
values. 
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APPENDIX D 


'I 


APPROXIMATE MIDCOURSE VELOCITY REQUIREMENTS 

The midcourse Av required to correct for injection errors is approximately equal 
to the uncorrected miss at the moon divided by the time to impact after the correction is 
made 


(Dl) 

m 

where M is the miss at the incoming asymptote to the moon. Equation (Dl) is valid be- 
cause the trajectory between midcourse and the moon is nearly a straight line (all gravity 
forces are small in this region). 

A correction time of 20 hours after injection is assumed, which gives T m a value 
of approximately 46 hours. The following analysis will determine the miss associated 
with the various injection errors. 


Error in Injection True Anomaly 

An error in true anomaly is equivalent to a rotation of the translunar ellipse about h 
by the error angle, so that 

M = rjAn | (D2) 

where r m is the target distance (Earth-to-moon radius) and Arj is the error in true 
anomaly. 


Error in Perigee Radius 

The energy and perigee radius can be expressed in terms of the orbital elements 


r 


P 


P 

1 + e 


(D3) 


E = — (e 2 - 1) 
2p 


(D4) 
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and the differentials are 


dE = - (e 2 - l)dp + -^ de 

2p 2 P 


dr _ (1 + e)dp - p de 
p o 

(1 + e) 2 


At constant energy, equations (D5) and (D6) give 


dp = 2e dr 


de = 


_ (e - 1) 


dr. 


The target true anomaly can be obtained from 


C0ST 7 m =- 


-P-- 1 

r m 


(D5) 


(D6) 


(D7) 

(D8) 


(D9) 


so that the equivalent error in true anomaly is 


dT ?m = 


JL - lW 

v r m / 


e 2 sin v r m esin ^m 


dp 


(DIO) 


In terms of dr , 

r 


drj 


dr p r r m,.2 


m 


r m e smT ?m 


— — (e - 1) + e + 1 
. P 


■] 


(Dll) 


and the miss at the moon is 


M = 


2 . 
e sin t] 


' Ar E_r!H(e 2 -l) + e 2 + ll 

IT, Lp -* 


m 


(D12) 
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Error in Injection Velocity 


At constant perigee radius, equations (D5) and (D6) give 


2p dE 

de = — — 

p(e + 1) 


(D13) 


dp = 


2p 2 dE 
p(e + 1)" 


(D14) 


Using equation (D10) results in 


2p dE 

d77 

pe(e + l)sin rj 


m 


1-JL 

v r my 


r (1 + e) 
m 


Energy is related to velocity and radius by 


(D15) 


r 2 
2 r 


E = - — 


(D16) 


At constant radius, 


dE = v dv 


(D17) 


and the miss at the moon is 


-2p r m v Av 

( ■ t) 

P \ m /_ 


pe(e + l)sin T? m 

r m (l + e) 6 



- 



(D18) 
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Error in Yaw Velocity 




f 


Yaw velocity is defined as 

v h =(?xf). R m = v e (h . H m ) (D19) 

An error in yaw velocity implies that the injection plane does not contain the target, and 
the corresponding miss is given by 


v e 


nv 


(D20) 


Numerical Sensitivities 

The following numerical values are assumed for the calculation of the numerical sen- 
sitivity coefficients: 


Time to impact, T m , sec 1.6x10 

Target distance, r m , ft 1.3xl0 9 

Eccentricity, e 0. 986 

Target true anomaly, rj , deg 170 

m 7 

Semilatus rectum, p, ft 4.26x10 

Velocity, v, ft/sec 3. 6xl0 4 

Central body force constant, p, ft^/sec^ 1.41xl0 16 

Circumferential velocity, Vg, ft/sec 3.6xl0 4 


These values and equations (Dl), (D2), (D12), (D18), and (D19), are used to obtain the 
following sensitivity coefficients: 


Injection error 

Midcourse Av 

True anomaly 
Perigee radius 
Velocity 
Yaw velocity 

137 (ft/sec)/deg 
0. 241 (ft/sec)/n. mi. 
5.21 (ft/sec)/(ft/sec) 
0.218 (ft/sec)/(ft/sec) 
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